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NET RADIATION METHOD FOR ENCLOSURE SYSTEMS 
INVOLVING PARTIALLY TRANSPARENT WALLS 
by Robert Siegel 
Lewis Research Center 

SUMMARY 

The net radiation method is developed for analyzing radiation heat transfer in enclo- 
sure systems involving partially transparent walls. One such system is an enclosure 
with windows in it. The conventional net radiation method has been well developed in 
the literature for enclosures with opaque walls. When one or more of the enclosure 
surfaces are transparent or partially transparent, radiation can directly enter or leave 
the enclosure. The enclosure equations are developed here to account for these effects. 
Specifically treated are gray and semigray enclosures. Another system of interest is a 
series of parallel partially transparent layers, such as used in flat-plate solar collec- 
tors. The characteristics of such a system are obtained here by the net radiation 
method which appears to be more convenient and systematic than the ray tracing tech- 
nique that has been used in the past. Overall transmission, absorption, and reflection 
relations are developed for windows consisting of any number of parallel layers with or 
without spaces between them, and that have differing absorption coefficients and differ- 
ing surface reflectivities. Also considered are the characteristics of combinations of 
parallel transparent layers and opaque walls. 


I PRODUCT ION 

The net radiation method has been well developed in the literature, such as in the 
textbooks (refs. 1 to 3), for calculating radiation heat transfer in enclosures. The- walls 
of these enclosures have been opaque. The object of this report is to show how the net 
radiation method can be modified and applied for systems involving walls that are par- 
tially* transparen t for radiant energy. Such systems would include enclosures with one 

*The term ’’partially" transparent means that there can be absorption of some of 
the radiation traveling within the wall with the rest being transmitted through the entire 
thickness. 



or more windows, or devices like flat plate solar collectors that consist of a series of 
parallel partially transparent plates with an opaque absorber surface behind them. 

The transmission characteristics of windows in the form of single or multiple glass 
plates have been derived by Whillier (ref. 4) using ray tracing techniques. It will be 
shown here how these results can be developed by use of net radiation concepts that 
provide a convenient systematic approach. 


SYMBOLS 


A 

a 

F 

F 

m 

N 

q 

r 

T 

a 


0-X c T 


surface area 

absorptance for path within layer 
radiation configuration factor 

fractional emission of blackbody in wavelength interval 0 to A c 
number of opaque surfaces in enclosure 
total number of surface areas in enclosure 
heat flux 

11 . iu 

surface reflectivity r at m surface of n layer 
absolute temperature 

overall absorptance of layer or system of layers 
a at plate a ks°rptance at plate in presence of transparent plate system 
/3 angle measured away from normal of surface 

e emittance of partially transparent layer; emissivity of opaque surface 

X wavelength of radiation 

A c cut-off wavelength for partially transparent material 

p overall reflectance of transparent plate system 

a Stefan- Boltzmann constant 

r • overall transmittance of transparent plate system 
Subscripts: 

e externally incident 

k k“ surface area 

i incoming radiation 
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I loss to outside 

0 outgoing radiation 

s for entire system 

Superscripts: 

1 for radiation incident on inside surface 
l long wavelength radiation A > A c 

(n) for system of n layers 

o for radiation incident on outside surface 

s short wavelength radiation X < A 

V % 

ANALYSIS 

Radiative Characteristics of a Single Layer 

To illustrate the transmission, absorption, and reflection characteristics of a win- 
dow for incident radiation, consider a single layer as shown in figure 1. The reflectivity 



Figure 1. - Net radiation quantities for ana- 
lyzing overall performance of partially 
transparent single layer. 


at each of the interfaces is r and the absorptance of the layer is a. The relations be- 
tween the outgoing and incoming energy fluxes at each interface can be expressed in 
terms of the reflection at the interface and the transmission across the interface. Using 
the net radiation concepts, the outgoing flux is written for each side of each interface to 
give 
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Vl =rq i,l + (1 ' r)q i,2 < la > 

V2 = (1 ’ r)q i,l + rq i,2 < lb > 

q o,3 =rq i,3 + ( 1 - r)q i,4 ( lc ) 

q o,4 = ( 1 " r ) q i,3 + rq i,4 < ld > 

The internal q.’s and q^s are related by using the absorptance of the layer to give 

q i, 2 = ^ " a ^ q o, 3 ( 2a ) 

q i, 3 = ^ " a ^%, 2 ( 2b ) ' 


The overall radiative behavior can be studied by considering energy incident on only one 

side of the layer, and superposition used when radiation is incident on both sides. Since 

it is the fractions of incident radiation transmitted, absorbed, and reflected, that are 

desired, the incident flux can be set equal to unity without any loss of generality. Then 

any q in the calculation will be a fraction of the incident flux. Hence, using q. 1 = 1 

and q. . = 0 along with equations (2), the q.'s are eliminated from equations (1) to 
1 ? ^ 1 

yield 


q Qjl =r + (1 - r)(l - a)q o 3 

(3a) 

q o,2 =1 ' r + r ( 1 - a ) q o,3 

(3b) 

*^ 0 , 3 = r ( 1 " a ) q o, 2 

(3c) 

q Q 4 = (1 - r)(l - a)q Q 2 

(3d) 


These are solved simultaneously to yield the reflectance p and transmittance r of the 
entire layer as 


p = V 1 = 


r 1 + (1 - 2r)(l - a)^j 


1 - r 2 (l - a) 2 


T = q o,4 


(1 - r) 2 (l - a) 
1 - r 2 (l - a) 2 


(4a) 


(4b) 
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The absorptance is then obtained from the overall heat balance 


a = 1 - p - t = ■ (4c) 

1 - r(l - a) 

Equations (4) can be directly applied when the incident energy is from a single 
direction. The r would be evaluated for that incidence angle and the a would be based 
on the path length traversed within the medium (allow for change in direction as radia- 
tion crosses interface of layer). If the incident flux is diffuse, the average p is obtain- 
ed by averaging over all angles as (ref. 1) 

fir/2 . 

p = 2 / p(/3) cos /3 sin /3 d/3 (5) 

J p=o 


where /3 is the incidence angle measured away from the surface normal. The p(/3) 
is given by equation (4a) with the r and a inserted as a function of /3. The r from 
equation (4b) can be averaged in the same manner and the average a then found as 

1 - p - T. 

In the enclosure theory that follows, the p and r for a layer refer to expressions 
such as equations (4a) and (4b), that is, relations that provide the overall performance 
of a window for incident energy. If the window is hot enough to radiate appreciably, the 
emittance of the window is also needed in the enclosure heat balances. According to 
Kirchhoff’s law, the e can be found by evaluating a for incident radiation having the 
same spectral distribution as blackbody emission at the temperature of the window. 

This assumes that the window configuration is such that a uniform temperature can be 
assigned to its entire volume. 


Enclosure Theory Equations for Gray Surfaces and Windows 

The net radiation method in the literature will now be generalized to include enclo- 
sure boundaries that are partially transparent to radiation. The usual restrictions of 

the enclosure theory apply which include diffuse emission and reflection, uniform fluxes 

» 

over each surface area and a uniform temperature over each area. A window such as a 
smooth glass plate will reflect primarily in a specular fashion so that the assumption of 
diffuse reflection is violated. However, within an enclosure there are extensive multi- 
ple reflections and the directionality of each reflection loses its importance in contribut- 
ing to the heat fluxes on the boundaries. Hence, the assumption of diffuse reflections is 
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usually satisfactory. The enclosure heat balances include possible emission from each 
window. The formulation here is only general enough to account for a window that can 
be characterized by a single temperature throughout its thickness. The formulation 
would require extension to account for a multilayered window with each layer at a differ- 
ent temperature. No extension is necessary, however, if the multilayered windows are 
kept cool by convection for example, so that emission from the layers is small relative 
to the transmission through them. 

A typical enclosure containing one or more partially transparent windows is shown 

i-L 

in figure 2. In general, let there be a total of N surface areas with the k area being 
a typical area. There are m opaque surfaces and N - m partially transparent sur- 
faces. 



Figured - Enclosure with N surface areas, some 
being partially transparent 


Outside 



Inside 

Figure 3. - Typrcaf nonsym metric window. 


Each window can be nonsymmetric about its center plane as shown in figure 3 . The 
superscripts o and i refer to the outside and inside surfaces. The a 0 and p° are 
the absorbed and reflected portions of the radiation incident from the outside. The e° 
is the emittance toward the outside and the t° is the fraction of the outside radiation 
that is' transmitted all the way through the window. A similar notation is used for the 
quantities at the inside surface (it will be found later that, for composite windows, 

T° = T*). 

For the m surface areas that are opaque the conventional net radiation equations 
can be written as: 

Overall heat balance: 
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^k “ q i,k 


(1 < k < m) 


(6a) 


Outgoing energy consists of emitted and reflected portions : 

%,k ' ‘ r k‘ rt[ k + (1 ' e k>«i,k (izksm) 
Incoming energy in terms of outgoing fluxes and view factors: 


N 


Vi,k=£%,i A i F j-k 

3=1 


(1 ^ k s m) 



(1 ^ k < m) 


(6b) 


.. (6c) 


These equations can be combined in various ways to provide other interrelations 
between the variables that are sometimes useful. Combine equations (6a) and (6c) to 
relate q and q 0 ’s, 


N 

(l=k = m) (7a) 

i=i 

Combine equations (6a) and (6b) to relate q, T and q 0 

(lsksm > ™ 

" k 

Combine equations (6b) and (6c) to relate the q 0 's and T, 

N 

%,k~ ^- e k^%,j F k-j= e k aT k (1 ~ ^ — m) « (7c) 

3=1 

Now consider a partially transparent surface. Using the quantities in figure 2 an 
overall heat balance yields 


% s V q i,k +q l,k' q e,k 


(8a) 
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The q^ could be supplied for example by heating wires within the semitransparent 
material or by convection on the outside surface. For cooling the is negative. The 
radiative flux leaving the inside surface of the window consists of emitted energy, re- 
flected incoming energy, and transmitted externally incident flux, 


<o, k “ e k ,7T k + fJ k q i, k 


T k q e,k 


Similarly the flux lost from the outside surface of the window is given by 


(8b) 


^,k = £ k aT k + p k'!e,k + ’Xk (8c > 

The incoming flux q. . is given by the same relation as equation (6c) 

1, K 

q i,k=Evi F k-i (8d > 

j=l 

Equations (8) are valid for m + 1 < k < N. 

Equations (8a) and (8c) are combined to eliminate q^ to give 

q k " e k crT k + < *e,k( 1 ~ p k) = %,k~ ~ 4:) ^ 

Some additional relations are found by combining some of the previous equations. Com- 
bine equations (8d) and (9) to yield 


N 

q k- £ k' TT k + 'le,k( 1 - p k)'%,k-( 1 -’k)£%,j F k-i (m + lsk = N) (10) 

j=l 


Combine equations (8b) and (9) to yield 


■ q k ““T |[ e k p k + ( X " T k) € k] CTT k " e k q o,k " [ p k( X " p k) " T k( 1 ~ 7 k)] q e,k| ^ 


where the relation 1 - -r^. - = o?^ = has been used for gray surfaces. Combine 

equations (8b) and (8d) to yield 
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( 12 ) 


N 


%,k " %,j F k-j 


3=1 


" e k crT k + T k q e,k 


The previous relations relate the q^, T^, q g and q o ^ for each of the enclo- 
sure surfaces. By writing these equations for each surface there are sufficient equa- 
tions to solve for the unknown values. For an enclosure with opaque walls this is dis- 
cussed in detail in references 1 to 3. One convenient approach is to eliminate the 
q^^ k's thereby obtaining a set of equations relating the surface temperatures T^'s, the 
surface heat flux additions q^'s, and the externally incident heat fluxes q e 's. 

As a first step the ^ is obtained as a function of and from equa- 
tions (7b) and (11) 



1 - e. 


% 


(1 < k < m) 


(13a) 



{[❖L + (‘-4)4K-4%-[4(i- 



(m + 1 ^ k < N) 
(13b) 


Equations (13) are then inserted into equation (7a) and into equation (10) to eliminate the 

N 

q^ . 's. The summation £ has to be made in two parts and equations (13a) and (13b) 
^k j= i 

used in their appropriate ranges. After simplification this yields for the opaque sur- 
faces, 
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]= m+1 

(1 ^ k < m) (14a) 


and for the partially transparent surfaces; 



Writing equation (14) for each k yields the final set of N simultaneous equations relat- 
ing the N surface temperatures, the N surface heat flux additions, and the N - m 
external fluxes on the windows. If 2N - m of these quantities are specified, the equa- 
tions can be solved for the remaining unknowns . 

The enclosure theory formulation contains within it the transmittance, reflectance, 
and emittance of each window. If a window is kept cool by convection so it does not emit 
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appreciably, then the window emission is not of significance in the analysis. The window 
can then be multilayered as the overall transmittance and reflectance can be obtained 
for such a system as will be shown in subsequent sections. If a window is hot as a re- 
sult of absorbing radiation or being convectively or internally heated and is emitting 
significantly, then the emission term involves the window temperature which is specified 
as uniform. This would restrict the formulation to simple windows for which a uniform 
temperature can be assigned. The net radiation method can be applied for more com- 
plex situations than those discussed here, including enclosures with multilayered win- 
dows, with each layer emitting and at a different temperature . 


Semigray (Two Band) Method for Enclosure Analysis 

The transmission of most window materials is very wavelength dependent having 
high transmittance in the visible and a portion of the near infrared and low transmittance 
at longer wavelengths. There is often a rather sharp transition between the high and 
low transmission regions. The cut-off wavelength at which this transition occurs can be 
typically at about 2.5 micrometers for ordinary window glass or 25 micrometers for a 
sodium chloride crystal. 

In the semigray method (ref. 5) or two-band approximation the property variations 
are accounted for by using two sets cif property values, one uniform value for wave- 
lengths shorter than the cut-off and a . second value for wavelengths longer than the cut- 
off. The notation e s , p s , and so forth, will be used for the short wavelength region 
0 < X < X c and 6^, p^, and so forth, will be used for the long wavelength region 
X < X < <*> . A common application of the semigray method is in dealing with devices 
exposed to solar energy. The incident solar energy is predominantly at short wave- 
lengths characteristic of visible radiation, while the reradiation is at long wavelengths 
characteristic of the device surface temperature . In the present discussion it will be 

assumed that the externally incident fluxes and their spectral distributions are known; 
s l 

hence, q e and are known. 

If all the surface temperatures are known, then equations (14a) and (14b) can be 
solved after small modifications. Two equations are written for each surface, one for 

4 

X < X and the second for X > X . In each of these equations the crT is replaced by 

L L 

the blackbody emission at T within the respective wavelength intervals 0 - X„ and 

4 4 c 

X c - 00 . These quantities are given by crT Fg_^ rp and oT (1 - Fq T ) where Fq-a T 

is the fractional emission of a blackbody in the range 0 - X c for a blackbody at T; the 
function Fq depends only on the product XT. 

Since all quantities on the right of each equation are given, this yields 2N equa- 
tions that can be solved directly for the unknown N values of q^ and N values of q£. 
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Then the heat flux being supplied to each of the k surfaces or windows by means other 
than radiation is 




( 15 ) 


The set of equations for two wavelength bands is given in the appendix. 

A complication arises in the simultaneous solution when it is the equilibrium tem- 
perature of a surface that is to be calculated when the q. for that surface is specified. 

s Z ^ 

In that instance the q k and q k are not known individually, and it is only known that 

their sum is q. . This introduces an additional unknown in the simultaneous solution 

K si 

since the T k must be found with the separate values of q k and q k still being involved 

in the equations as unknowns as it is only their sum that is known. The equations given 

in the appendix still apply. There are equations (14a) and (14b) written for each X 

range which gives 2N equations, and equation (15) written for each of the N surfaces. 

In general, let there be p surfaces that have q k specified; hence N - p surfaces 

have T k specified. The total of 3N simultaneous equations can then be solved for the 

following unknowns : p values of T k , N-p values of q k , N values of q k and N 

values of q k . These equations contain, however, the term ctT^Fq_^ rp and the tern- 

0 

perature is not contained in Fq_^ rp in a simple way. Hence, an iterative procedure 

is required to determine the unknown T values. 

If the property variations with wavelength are such that many more than two wave- 
length regions are required to adequately account for the property variations, then the 
net radiation equations can be written in spectral form and the fluxes then integrated 
over wavelength, or a summation made over a finite number of wavelength bands. This 
is an extension of the two band procedure and has been described in textbooks such as 
reference 1 (chapter 10); hence, the procedure will not be further detailed here. 


Transmission Relations for Various Types of Windows 

Single layer with different surface reflectivities . - A geometry that provides a use- 
ful building block to obtain results for more complex windows is a single layer with a 
different reflectivity at each surface. A dual subscript notation will be used for the 
surface reflectivities: For a group of several parallel layers, the first subscript desig- 
nates the layer number and the second designates the upper or lower surface of that 
layer (the layers will be drawn horizontally). Thus, r k 2 would refer to the second or 
lower surface of the k^ 1 layer. Referring to figure 4 the net radiation equations for 
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Figure 4. - Single absorbing layer with 
different reflectivity at each surface. 


the first of several possible layers are written as 


Vl =r ll q i,l +(1 ‘ r H )< ‘i,2 < 16a > 

<Jo,2 = f 1 * r ll)«i, 1 + ' < 16b > 

%,3= r 12 q l,3 + (1 - r 12> < li,4 < 16c > 

%,4 = < 1 - r 12> q i,3 + r 12% > 4 (16d > 

The absorptance of the medium is used to relate the internal q^s and q Q , s to yield 

< M = (1 ' a l )q o, 3. ' (l7a) 

q i,3 = t 1 " a l^ q o, 2 ( 17b ) 


To consider the interaction of the layer with radiation incident on the upper surface, 

let q. 1 = 1 and q. A = 0. Using these relations and equations (17) to eliminate the 
1 ; 1 * 
q/s from equations (16) give 


q o,l =r ll + < 1 - r ll>< 1 - a l>V3 

(18a) 

lo,! 5 * 1 ' r ll> + r ll( 1 ‘ a l> q o,3 

(18b) 

V3 = r 12 (I ' a l )q o,2 

(18c) 

% j 4 = (l- r 1 2Kl- a i)q 0 ,2 

(18d) 
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The solution yields the overall transmittance and reflectance for radiation incident on 
the upper surface as 


T 11 “ %, 4 ~ 


(1 - r 11 )(l - r 12 )(l - a x ) 
1 - r ll r 12< 1 - a l> 2 


P 11 =%,1 


r ll - | ' r 12< 1 ~ 2r ll> (1 ' a l>‘ 
1 - r ll r 12< 1 - a l> 2 


(19a) 


(19b) 


The subscript 11 indicates that the and are for radiation incident on the first 
surface of the first of several possible layers. The overall absorptance for the layer is 
obtained from the energy balance 


O-ll = ! - PH - Tu 


Substituting the previous expressions for and p^, this simplifies to 


a 
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_ a l^ " r ll^[ 1 + r 12^ ' a l- > ] 
1 - r u r 12 (l - a,) 2 


(19c) 


If radiation is now considered incident from below (that is on surface 12), the only 
change in the analysis is that r^ and r ^ are interchanged. The overall transmit- 
tance, reflectance, and absorptance are designated by Tj 2 , P^, a-j 2 an< * * he y are 
given by 


t 12 


(1 “ 

1 " r 12 r ll (1 " a l) 2 


'12 


a 


12 


r 12 + r ll( 1 - 2r 12 ) ^ 1 ~ a l } 

1 - r 12 r ll( 1 - a l) 2 
_ a lU - r 12 )[l + r n (l - a t )] 
1 " r 12 r ll( 1 - a l) 2 


(20a) 


(20b) 


(20c) 
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This yields the relation = t-^ so that the second subscript for r can be omitted 
and the overall transmittance of the layer called Tj. The emissivity of the layer can 
be found according to Kirchhoff 's law by evaluating a using window properties for in- 
cident radiation that has the same spectral distribution as blackbody radiation at the 
window temperature. 

Two differing layers with all reflectivities different . - The behavior for a single 
layer can be used to obtain the performance for two layers as shown in figure 5. The 



Figure 5. - Two differing absorbing layers with 
all surface reflectivities different. 


q 0 's are written in terms of the q^s using the overall characteristics of single layers, 


V 11 * «1, ll p ll + % 12 t 1 < 2la > 

V 12 =“i, 12012 + Vl T l (21b) 

%,21= t ll,21 p 21 +< )l,22 T 2 < 21c > 

V 22 " *11,22022 + 9 i,21 T 2 < 21d > 


To obtain the overall characteristics of the two layer system for radiation incident from 
above, let q. = 1 and q^ 22 = 0* Also, from figure 5 

q i, 12 = %, 21 
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and 


q i,21 - V 12 

These relations are used to eliminate the q^s from equations (21) to yield 


%, 11 =p ll + q o,21 T l 

(22a) 

%,12 =q o,21 p 12 + T 1 

(22b) 

%,21 = %,12 p 21 

(22c) 

%, 22 = q o, 12 t 2 

(22d) 


Now introduce the notation r g j, r g 2 and p j, p^. The subscript s refers to 
the entire system . Hence, p g j and p g 2 are, respectively, the reflectances for radia- 
tion incident on the first and second surfaces of the entire system. The simultaneous 
equations (22) are solved to yield 


Ti t , 


T sl " ^,22 “ — 


1 2 

p 12 p 21 


(23a) 


2 ^21 

p sl =q o, 11 =p ll + T 1 Y 


(23b) 


P 12 P 21 


“si " 1 ‘ T sl ‘ p sl " 1 " P 11 


1 - 


p 12 p 21 


(t 2 + r l p 21 ) 


(23c) 


By changing subscripts to reverse the order in which the radiation passes through 
the layer surfaces, the overall system characteristics are obtained for radiation inci- 
dent from below on the two layer system, 


t s2 


t 2 t 1 

1 ‘ p 21 p 12 • 


(24a) 
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(24b) 


p s2 = p 22 


2 

T 2 P 12 

+ 

1 " p 21 p 12 


t 2 

“s2 = 1 “ p 22 “ ~ ( T l + t 2 p 12^ ( 24c ) 

1 - p 21 p 12 

A comparison of (23a) and (24a) reveals that r g j = T g 2 so that the system transmittance 
can be called simply r g . 

Two layers with no space between them . - Two adjacent layers with differing prop- 
erties are shown in figure 6. The reflectivity at the common interface between the two 


r n 

r l-2 


r 22 

Figure 6. - Two differing adjacent absorbing 
layers. 


| Layer 1 

3 1 

| Layer 2 

a 2 | 


media is called r^_ 2 » The configuration in figure 5 will correspond to this case if we 

let r ^2 = r i_2 anc * r 2l = ® (° r us i n S r 2 i = r i-2 anc * r 12 = 0 yield the same final 
result). The overall transmittance of the system can be computed by using equation (24a) 


t 1 t 2 

1 " p 21 p 12 


(25) 


where from equations (19) and (20), 

(1 - r n )(l - r 1-2 )(1 - a x ) 

J Tj = 

1 " r ll r l-2^ 1 " a l^ 
t 2 = ( J - r 22 )(1 ' a 2> 


(26a) 

(26b) 
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(26c) 


r l-2 + r ll< 1 - 2r l-2 )(1 ~ a l )2 
1 " r ll r l-2^ ~ a l> 2 

P 21 = r 22^ ” a 2^ ( 26d ) 


The expressions in equations (26) can also be used in equations (23b and c) and (24b 
and c) to compute p g p a s p p g and for the system of two immediately adjacent 

layers. 

A general system of N layers . - The overall characteristics of one layer were 
used to build up the overall behavior of a two layer system. This procedure can be 
continued and the behavior for one layer and two layers can be combined to yield the 
behavior for three layers; two layers and two layers can be combined to yield the be- 
havior of four layers; etc. To accomplish this it is desired to have relations for deter- 
mining the performance of a combined system of n layers and k layers, where 
n + k = N. 

Referring to figure 7 let T g n+ ^ be the system transmittance for the entire system 
and let be the transmittances for n and k plate systems, respectively. 

The quantity p^ is the reflectance for a system of n layers by itself for energy in- 
cident on the second (lower) surface of that entire n layer system. Similarly, p^ 


Layer 1 


a l 
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Figure 7. - System of N = n + k plates. 
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is the overall reflectance for a system of k layers by itself for energy incident on the 
first (upper) surface of that entire system. 

Comparing figures 7 and 5 indicates that the results for the two layer system can 
be applied if we use the overall properties of the separate n and k plate systems in 
figure 7 to replace, respectively, the properties of the first and second single plates 
in figure 5. Then from equations (23a), (23b), and (24b), 


. (n) (k) 

(n+k) _ T s T s 


1 . n(n) n (k) 


p s2 p sl 


(27a) 


2 i>0O 

(n+k) _ (n) + T (n T P sl 


p sl =p sl 


i - ofiM? 


(27b) 




2 >> 
1 p s2 


i - $$ 


(27c) 


IV -i. 

If the spacing between two plates is zero, say the m l and m + l 51 layers are in per- 
fect contact, then let r m2 = r m _ (m+1) and let r (m+1)1 = 0. 

As an example, compute t_ for a three layer system. Letting n = 2 and k = 1 
in equation (27a) gives 



t (2) t U) 
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where 


Then 


t ( 2). = T 1 T 2 
S 1 ‘ p 12 p 21 



“ p 22 


+ 


2 n ' 
t 2 p 12 


1 " p 21 p 12 



" p 31 


r(S)_. 


U " p i2 p 2lM 1 ' p 22 p 31 " ~ 


t 1 t 2 t 3 T 1 T 2 T 3 

T 2 p 12 p 3l\ ^ " p 12 p 21^ 1 " p 22 p 31^ “ r 2 p 12 p 31 


p 21 p 12> 


(28) 


where the required characteristics for single layers are given by equations (19) and (20) 


as 


=■ 

m 


a - r ml )(l - r m2 )(l - a m ) 




1 - r ml r m2 (1 - a m> 


ml 


1 ' r ml r m2< 1 - a m> 2 


V m = 1, 2, 3 


p m2 ' 


r m2 * r ml< 1 - 2r m2>< 1 - a m>‘ 

1 ' Vmll 1 - a m> 2 




As a special case if there is no absorption in any of the plates, a m = 0 and, if in addi- 
tion, all the surface, reflectivities are equal, then equation (28) reduces to 
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where 


T = 


(1 - r)' 


1 - r 


1 - r 
1 + r 


P = 


r + r(l - 2r) 


1 - 


2r 

1 + r 


If r and p are inserted into the expression for 



it simplifies to 


(3) 1 - r 

s 1 + 5r 


This agrees with the general result in reference 6 for a series of n layers with a space 
between each of them and with all surfaces having the same reflectivity r 


T (n) = 1 ~ r 

s 1 + (2n - l)r 


(30) 


Systems of Parallel Partially Transparent Layers and Opaque Plates 

A system of partially transparent layers with an absorbing surface below them . - 
This is a typical solar collector geometry and is illustrated in figure 8. It is desired to 
know how much of the incident solar radiation is absorbed at the surface N. For a unit 
flux incident on the first surface of the first layer, q i ^ = 1, the fraction absorbed at 
the plate can be obtained from 


^at plate ^i,Nl ” %, N1 

The net radiation method with q. ^ = 1 yields the following relations: 

n - n /-,( n ) , Jn) 

q o, (N-l)2 - % (N-l)2 p s2 + T s 


(31) 


(32a) 
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Figure 8. - System of n absorbing layers followed by absorber 
plate. 


where n = N - 1 layers, 
tion (31) yield 


q o,Nl = r Nl q i,Nl 

(32b) 

q i,Nl = V (N-l)2 

(32c) 

q i, (N-l)2 = %, N1 

(3 2d) 

Solving for q iN1 and % m 

and substituting into equa- 

_ (i - r N1 )4 n) 

at plate 

1 - *nipS 

(33) 


As a simplified case, if there are n nonabsorbing layers all having the same surface 
reflectivities r, then using equation (30) and the fact that p = 1 - t for a nonabsorbing 
system, yields, 


T ( n ) = 1 - l 

s 1 + (2n - l)r 




T (n) = 2nr 
s 1 + (2n - l)r 
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If these relations are substituted into equation (33), the fraction absorbed by the plate 
becomes 


a 


(1 “ r N1 )(l " r) 


(1 - r Nl )(1 " r) 


at plate 


1 + (2n - l)r - r N1 2nr (1 - r) + 2nr(l - r M1 ) 


( 34 ) 


Nl' 


where n = N - 1 layers. 

Some additional results can be obtained quite easily for the systems shown in fig- 
ure 9. For the situation in figure 9(a), = r 12 = r 21 from equations (19a) or (20a) 
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Figure 9. - Single absorbing layer inter- 
acting with opaque surface. 


and Pg*2 = from equation (20b). Making the substitution in equation (33) yields 


a 


(1 “ r Nl^ " r i2^0 ~ r u)0 " a l^ 


at plate o 2 

1 " r 12 r ll^ " a l^ " r Nl r 12 " r Nl r ll^ 1 " 2r 12lO " a l^ 


(35) 


H r u = r ^2 = r as a special case, equation (35) reduces to 


a 


(1 - r N1 )(l - rHl - a x ) 


at plate o o 9 

1 - r (1 - a - r N1 r - r N1 r(l - 2r)(l - a x r 


(36) 


If in addition the layer is nonabsorbing, a^ = 0, the result simplifies to 
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(37) 


“at plate 


(1 - r N1 )(l - r) 
1 + r - 2r N1 r 


For the situation in figure 9(b) the a ^ plate can be obtained from equation (35) by 
letting r 12 = r i_N and r Nl = This yi elds 


. » - rj.jjHl - r, 1 )(l - aj) 
“at plate 9 

1 " r l-N r ll (1 ‘ 


(38) 


Transparent layer between two opaque plates . - As a further illustration of the ap- 
plication of the net radiation technique, consider the effect on the heat exchange of plac- 
ing a transparent nonabsorbing layer between two parallel plates as shown in figure 10. 


q 



q 


Figure 10. - Effect of intermediate transparent 
layer on exchange between two opaque par- 
allel plates. 


The net radiation equations are written as follows: 

At plate 1 

«i,4 +< l = V4 < 39a > 

V4 = £ l ffT l + (1 ' £ l )q i,4 < 39b > 


24 




At plate 2 


For the transparent plate 


+ v i 

%, ! =«2» t 2 +(1 ‘ £ 2)«i,l 


(40a) 

(40b) 


q o,2 = q i,3 T+q i, 2 P 
q o,3 =q i,2 T+q i,3 p 


(41a) 

(41b) 


where p and t are the overall reflectance and transmittance of the transparent plate. 
The connection between the q^'s and q^s is given by 


q i, 3 ~ %, 4 
q i, 1 = %, 2 


q i,2 = q o,ll 


q i,4 = q o,3j 


(42) 


The relations in equations (42) are used to eliminate the q/s from equations (39) to (41) 
to yield 


%,3 + « = %,4 
%,4 s£ l‘’ T l +(1 - £ lK>,3 
1o,2 = 1 + Vl 
%, 1 = € 2 cT 2 + < 1 ' ' 2 * 10,2 
%,2=1o,4 T+ 1o,l p 
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The 2 and q c 3 are eliminated to yield 

q 0; 4=€iCTj t (l-£ 1 )(q 0j4 -q) 

Vl = £ 2 aT 2 + 0 - € 2>('1 + %,1> 

1 + %,l=9o,4 T+ Vl P 

Vr‘i^,i Ti V# 

These relations are solved for q which is the heat flux being transferred from the low- 
er plate maintained at Tj to the upper plate at Tg. The result is simplified by using 
the relation p = 1 - r valid for a nonabsorbing layer to yield 



For a single transparent layer with surface reflectivities r, equation (4b) yields 

1 + r 

so that by substituting into equation (43) the final result is obtained that 



€ 1 e 2 1 ‘ r 

The 2r/(l - r) term is the contribution provided by the transparent plate. 


CONCLUSIONS 

The net radiation method has been generalized to situations where partially trans- 
parent layers are present, such as an enclosure having one or more windows. The 
utility of the method is further illustrated by applying it to determine the overall trans- 
mission and reflection characteristics of systems of parallel plates with or without par- 


26 



allel opaque plates. These configurations arise in solar collectors, solar stills, and in 
coated surfaces where the coatings are thick relative to the wavelength of the radiation 
so that interference effects do not have to be included. The method is found convenient 
to apply and provides a systematic approach that is not difficult even for complex sys- 
tems. Most partially transparent plates have spectrally dependent properties. The net 
radiation equations can be applied for various wavelength. regions and the energy quanti- 
ties then summed over wavelength. The semigray (two band) approximation has been 
specifically discussed here which utilizes two spectral regions in each of which the 
properties are constant. The utilization of multiple wavelength bands proceeds in the 
same fashion as has been well documented in the literature for systems with opaque 
walls . 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 17, 1973, 

770-18. 
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APPENDIX - ENCLOSURE EQUATIONS FOR TWO BAND APPROXIMATION 


The equations for the short wavelength region are for the opaque surfaces, 
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and for the partially transparent surfaces, 
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The equations for the long wavelength region are for the opaque surfaces, 
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and for the partially transparent surfaces, 
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The total heat fluxes are found by summing the fluxes for the two wavelength regions, 


q k = ^ + q k (1 — k — N) 


(A5) 
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